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Abstract 

Simple equations for the linear response of layered superconductors with d- 
wave symmetry of the order parameter are derived by means of kinetic equa- 
tions for Green's functions. Responses to solenoidal and potential electric 
fields have different frequency dependencies. The damping of plasma oscilla- 
tions of superconducting electrons is determined by dielectric relaxation and 
is small. Relaxation of branch imbalance determined by elastic scattering is 
large enough to make the Carlson-Goldman mode in d-wave superconductors 
overdamped. 
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Many evidences for (i-wave symmetry (or near <i-wave symmetry) of the superconducting 
order parameter in layered high-T c superconductors were given last years by measurements 
of the Josephson effect JI|-£§, by microwave experiments j4|, by high-resolution angle- resolved 
spectroscopy || and by other methods. On the other hand, theoretical studies confirm a 
compatibility of many experimental data with <i-wave symmetry of the superconducting gap 
in high-T c superconductors (see |],[7| and references therein) . Many properties of the <i-wave 
superconductors are expected to be different ^from those of conventional superconductors, 
especially the effects related to quasiparticles and their relaxation, since due to the nodes of 
the (i-wave order parameter the quasiparticle density is never exponentially small. Here, we 
study theoretically the effects related to electric field, and collective oscillations in <i-wave 
sup er conduct ors . 

Typically, calculations of the linear response of the superconductors assume the response 
to a solenoidal (transverse) electric field which can be expressed in terms of time deriva- 
tive of the vector potential, paying less attention to a potential (longitudinal) electric field 
which is determined as the gradient of the scalar potential. However, the linear response of 
superconductors depends on the origin of the electric field. Electric field created by varia- 
tions of the current density and magnetic field in time, and related to the Faraday's law, 
creates in a superconductor only the perturbations of the electronic distribution with the 
antisymmetric angle dependence. Such a perturbation is limited by the momentum relax- 
ation, like in a normal metal. On the other hand, the potential electric field which is related 
to perturbations of the charge density and to the Coulomb's law creates, in addition, the 
branch imbalance p|,p|JTT|, i. e. the difference between the densities of electron-like and 
hole-like quasiparticles. Thus, relaxation of the branch imbalance is involved in the linear 
response of a superconductor to the potential electric field as well. The branch imbalance is 
characterized by the gauge invariant scalar potential fi = (l/2)(dx/dt) + $ where x is the 
phase of the order parameter and $ is the electric potential. Potential \i can be interpreted 
as the nonequilibrium shift of the chemical potential of the normal carreers. Furthermore, 
the superconducting momentum, P s = (1/2) Vx — A where A is the vector potential, plays 
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a role of the gauge invariant vector potential (see [II] and references therein). Then the 



electric field is expressed in terms of the gauge invariant potentials as 

E=-V„ + ^. (!) 

So, two contributions to the electric field in ([!]) produce different kinds of the perturbations 
of the electronic distribution in a superconductor. We shall consider the first term of this 
expression as the potential part of the electric field. The second term related to the variations 
of the superconducting current in time we shall call the solenoidal part of the electric field. 
Strictly speaking, the second term in ([J) does not satisfy the usual definition of the solenoidal 
field because its divergence does not vanish exactly in all the cases. Nevertheless, we call 
it solenoidal or transverse field because it is needed to describe the purely solenoidal field, 
and it is related to the response to the electric field created by the time dependent magnetic 
field. 

The potential electric field must be taken into account in the problem of the linear 
response since it is important in collective oscillations and it appears in nonuniform and 
anisotropic systems, even if externally applied field is purely solenoidal. In order to study 
the problem we calculate a linear response of <i-wave superconductors to the gauge invariant 
vector and scalar electromagnetic potentials, and derive simple and physically transparent 
expressions for the charge and current densities from the equations for quasiclassical Green's 
functions using the nonequilibrium approach by Keldysh [K| . Such an approach enables us 



to take into account both momentum and branch imbalance relaxation rates. The latter 
enters the generalized conductivity describing the response to the potential electric field. 
Then we use these expressions to study plasma oscillations of superconducting electrons and 
the Carlson-Goldman mode in <i-wave superconductors. 

To calculate the linear response of layered superconductors with o?-wave pairing, we start 
with the equations for Green's functions in Keldysh technique, using the slightly modified 
approach by Larkin and Ovchinnikov [|TjJ. We use two different ways to describe layered 
superconductors. In the first approach we use a continuous representation, considering an 



anisotropic metal with the band motion of electrons in the direction perpendicular to the 
layers such that t± ^> u, where t± is the overlap integral describing the electron spectrum in 
perpendicular direction, e± = 2t± cos dp±. Here d is the lattice constant in the perpendicular 
direction, and v the momentum scattering rate along the layers. In the second approach we 
use similar equations [TB||H| | for layered superconductors in the discrete Wannier represen- 
tation, considering the hopping conductivity regime between the layers, i. e. t± <C v, which 
corresponds to the case of Josephson interlayer coupling. The second approach bears some 
similarity to the interlayer diffusion model |l5j in which the interlayer coupling is mediated 
through incoherent hopping processes with t± neglected. In both cases we assume a ci-wave 
superconducting order parameter: thus we do not address the question of the microscopic 
nature of the interaction resulting in such a symmetry. 



To derive the equations in the continuous representation we subtract, similar to ||12|| , from 
the equation for matrix Green's function in Keldysh representation its conjugated equation. 
Then we integrate the resulting equation over £ = p|/2m — sp, where py is the component of 
the momentum parallel to the layers. Thus we obtain equations for the retarded (advanced) 
Green's functions, g R ^ A \ and for g K , which is related to the electron distribution function. 
Each of these functions is a matrix in Nambu space and depends on coordinates, energies, 
perpendicular component p± of the momentum, and on the angle of pn. 

In the linear approximation with respect to the external perturbation, the equa- 
tion for the anomalous Green's function, g( a \ defined by g K = g R (s, e') tanh (e'/2T) — 
g A (e } e') tanh (e/2T) + g^(e, e') has the form 

vV ff (a) - [e + a z + A(<f>)i<T y ]gto + g {a) [e^a z + A{<f>)ia y ] - 

( E R g (a) _ gR^a) + ^a) g A _ = ^[(vP^^ + figA _ g R^V s a z + //)]. (2) 

Here v is the electron velocity at the Fermi surface, A(0) is the amplitude of the order 
parameter, a = tanhe + /2T — tanhe_/2T, Pauli matrices, and the unperturbed 

retarded and advanced Green's functions, g R ^ A \ depend on shifted energies e + = e + uj/2 
and e_ = £ — uj/2, respectively. The self-energy parts, S 1 , are given by 



J-n/d An/a Jo Zn 

where i =R, A, or K; v is the elastic scattering rate in the normal state. Strictly speaking, 
(§) describes the impurity scattering in Born approximation, but it can be applied also to 
elastic scattering by phonons, since the related self-energy part acquires the form (|3|) when 
one neglects phonon frequencies in comparison with electron energies in delta-functions 
describing the energy conservation law in scattering processes. Using Born approximation we 
neglect low-energy quasiparticle bound states created by impurities (see |16|] and references 



therein), and, hence, our results are applicable provided typical energies of quasiparticles 
are larger than the bandwidth of the impurity induced bound states, T > \J H\v. 

Now we consider the momentum dependence of the scattering rate v. As it will be seen 
below, the in-plane scattering results in pair-breaking similar to magnetic impurities in s- 
wave superconductors, while the interlayer scattering does not affect the gap. In addition, 
from the conductivity anisotropy data in high-T c superconductors one may expect, that 
corresponding components of v have different temperature dependencies. Having this in 
mind we consider a simple model for momentum dependence of u, which takes into account 
different scattering rates in different directions: 

v(p±, <P;P±, 0') = Vi + VxH<t> ~ </>')• ( 4 ) 

Here z/j describes the isotropic scattering, and u± is related to the scattering in the perpen- 
dicular, interlayer, direction. 

Using (|3D and we obtain for the unperturbed retarded (advanced) Green's functions in 
(|) the implicit relations g R ^ = a z a R ^ + ia y b R{ - A \ where a R{ -^ = (e + iui(a R ^)^/2)/^ A \ 
foR(A) _ A(0)/ £ R ( A ) . The brackets (...) mean averaging over variables mentioned in the 



subscript, and = ±y (e + iv i {a R ^) ( f ) /2)' 2 — A(0) 2 . The equations for perturbations of 
gR(A) can OD ^ a i ne( j from (^) by replacing a by 1, changing all the superscripts for R(A), 
and omitting the two last terms in the collision integral. 

The main differences in the equations for the Green's functions in the discrete represen- 



tation (see [|T4]]) are the following. Green's functions become matrices in layer indices, and 



the potentials depend on the layer index as well. The first term in the l.h.s. of (|2]) in the 
discrete representation is to be replaced by t±J2i=±i (A n n+i9n+im ~ 9nm+i^m+im) , describing 
the interlayer interaction, with A nm = cos (%„ — Xm)/2 + io z sin (% n — % m )/2. Averaging in 
the collision integral is performed over the angle (f> only. 

We solve the linearized equations for Green's functions for the case of smoothly varying 
perturbations |qv| <C Pi, i- e. when changes of all variables along the distance of the order 
of the mean free path along the layers are small. This case covers the most interesting range 
of frequencies, because characteristic values of 1/q are determined either by the magnetic 
penetration lengths (at low frequencies) which are typically larger, than the mean free path 
in high-T c superconductors, or (at high frequencies) by the skin-effect length, which is also 
large provided the frequency is below the range of the anomalous skin-effect. 

Perturbations of charge density are determined by Tr(g^) integrated over energies. 
Quasiparticle current densities in the directions parallel and perpendicular to the conducting 
layers are proportional to integrals of Trcr 2 (vg( a )) and of Tra z (v z g^) , respectively. Super- 
conducting currents are determined by similar terms with retarded and andvanced Green's 
functions in the relation between g K and g( a \ We calculate current and charge densities 
assuming the clean limit, T c 3> z/j, since in the opposite dirty limit a superconductor is in a 
gapless state. In the case of frequencies u much smaller than the amplitude of the gap, A 
the linear response can be presented in a simple and physically transparent form: 

K 2 

- iujp = -iuj>y—ij, + (a 2 iq 2 + cr 2t k 2 )iJ, + w^uqP, + a u kP t ), (5) 

47T 

c 2 

ii = -rv2 P i ~ i ( UJ Voi'Pi + ^ljq/i), (6) 

c 2 

jt = -r~v2 P t ~ i&VotPt + (Titkfi), (7) 

where P; and P t are parallel and perpendicular to the layers components of P s , q and 
k are parallel and perpendicular components of wave vector, k _1 is the Thomas-Fermi 
screening radius, and are the penetration lengths for a superconducting current parallel 
(perpendicular) to the layers. 
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The first terms in (H0) describe the supercurrents, while the last terms are related to 
the quasiparticles. Taking into account that the electric field expressed in terms of the 
gauge invariant potentials is given by (|I]) we see, that the simple expression j = <5~E for 
quasiparticle contributions to the currents is not valid: generalized conductivities a na are 
different for the contributions from scalar and vector potentials to the electric field. This 
implies different responses to the potential and to the solenoidal electric fields. Then note 
that according to (||) the changes of the charge density are determined by time variations 
of the potential /i, which is related to the electrone-hole imbalance (cf. P,|TT|), and by space 
variations of the quasiparticle currents. Equation (|5]) plays a role of the continuity equation 
for normal carreers. 

The factor 7 and the generalized conductivities a na {n = 0, 1, 2 and a = I, t) depend on 
frequency: 

7 = 1+/ de- — . , 8 

J -00 {00 + iVb) de 

1 r ^ ^(lel-lA^Dal-^ dn^ 

Vna = -VNa— / d£{— - w , ■ v_ "J— )«*J (9J 

r a J -00 (u + iv a )(u + ivb) de 
Here a^ a is the normal state conductivity in direction a, np is Fermi distribution function, 



and ao = ej \/e 2 — A 2 ; ^ = 1/^ = 1/^ and r t = l/(i/i + u±) are the momentum scattering 
times of electrons in the normal state for longitudinal and transverse directions, and v\ = 
^(oo)(/> an d i>t — V\ + v±/a are energy dependent effective scattering rates of quasiparticles 
for corresponding directions. Finally, i/& = ui(A 2 ((f))aQ/e 2 ) ( j ) is the effective branch imbalance 
relaxation rate. It is well-known that in s-wave superconductors the branch imbalance relaxes 
via inelastic scattering, spin-flip scattering or due to anisotropy of the order parameter (for 
a review see JTl| ) . In the case of <i-wave pairing the elastic scattering is a main source of the 
branch imbalance relaxation. 

In the studies of the linear response to a solenoidal field, the conductivity a^i in ([]) is 
usually calculated. Our result for cr « agrees with the Born limit of the general expressions 
for the conductivity obtained in 0. 

Note that the conductivities in the transverse direction are determined by contributions 



both from intralayer scattering and by interlayer scattering, so that the effective scattering 
rate for the conductivity in transverse direction is larger than the effective scattering rate 
for the in-plane quasiparticle current. 

Solution of the discrete equations in the limit of small phase differences between the 
neighbouring layers gives results similar to (|§-§5|); it can be obtained from (^-0) substituting 
Pt for (xn — Xn-i)/d, and v z for 2t±d with u± = and Vi = v in (||-|), so ^ na ^ conductivities 
in both directions are determined by the same scattering rate. 

Now we discuss the limits of low (A ^> T) and high (A <C T) temperatures. 

1. A ^> T: an important distinction from s-wave superconductors is that the conductiv- 
ities in (H) are not exponentially small. 

Consider, first, the linear response to the electromagnetic wave. For simplicity, in order to 
get explicit expressions we consider the simplest angular dependence of the gap parameter 
with the rf-wave symmetry A = A o cos20 . At low temperatures, (ao}^ = s/Aq, and 
the characteristic times for quasiparticles averaged over energies are f\ = r^(A /2T) « 
r b 3> Tj. The relative density of superconducting component in c?-wave superconductor is 
N s = l- (T/A ) In 4 = (A(0)/A(T)) 2 « l.The factor 7 is also very close to 1, 7 = 1 + 2iujT t 
for uo <C 1/tj, and 7 = 1 — (T/A ) In 4 for u 3> l/r^. This leads to a smaller contribution 
of scalar potential \x (though not exponentially small as in s-wave superconductors) , and we 
may omit the diffusion contribution to the quasiparticle current densities. Then the current 
densities can be presented as 
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J" = TTT^^ - iucr Na P a R a {u) (10) 

R [°° X< ^ X (11) 

Jo (x — iujfi) cosh 2 x ' 

Rt - - f dx t XdV t 5-, (12) 

7ir t Jo Jo [x — iujTi + v\Ti\J 1 — y 2 ) cosh x 

According to ( ]T0| - [TT] ) , at ujt\ <^ 1, i. e. when scattering is important, the decrease of the 
normal carreer density is compensated by the decrease of the scattering rate of quasiparticles 
in comparison with the normal state by the same factor oc T/A. Furthermore, Ri = 1, and 
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the quasiparticle conductivity along the layers even at low temperatures is the same as 
it would be in the normal state at this temperature. At higher frequencies (ufi ^> 1) 
Ri = i(l — N s )/(uTi), scattering is not important and the current density corresponds to a 
free motion of all electrons. 

Consider now the conductivity in the transverse direction. Note that due to the depairing 
action of the in-plane elastic scattering, these processes contribute to the conductivity in 
the transverse direction. If u± <C 1/fj or u± <C uj, we find Rt = -R;(r//r t ). Thus in this 
case the transverse conductivity is determined by the in-plane scattering. At u±_ 3> 1/tj and 
u± ^> uj, we obtain R t oc (T/A ) hiUj_/(u 2 + 1/f/ 2 ) is small and can be neglected. 

Using equations (|T0l - |i2l) one can easily calculate the surface impedance of a d-wave su- 
perconductor. For a surface parallel to the layers we obtain 



c = 



1 - (uj/uj ) 2 - iR t uJT t + (q\ t ) 2 ^ 



c ^ (1 - (uj/uj ) 2 - iR t uT t ){\ - iRiUJTi) ' 

where u = cj \ is the frequency of plasma oscillations for an electric field perpendicular to 
the layers. 

We conclude that, in spite of the large quasiparticle conductivity equal to the normal- 
state conductivity, the damping terms at low temperatures are always small, because the 
scattering of quasiparticles is important only at uj < 1/f; 1/t|. 

Now we discuss the spectrum of free oscillations, which can be calculated inserting (|^- 
0) into the Maxwell equations. The spectrum of the weakly damped plasma mode in the 



long- wavelength limit is given by an expression similar to that of the case of s-pairing ]14j : 

" = "H rrk? — lftwr, j (14) 

The last term in (13) describes damping. For small frequencies limit it is determined by 



the dielectric relaxation frequency ^on, which is rather large. Nevertheless, the plasma 
oscillations survive, because the damping is determined by only at frequencies uj < 
T/(rA) <C 1/t, and becomes small at uj ~ ujq > T/(rA). 
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2. A <^T : at high temperatures (but outside the gapless regime A < v{) the branch im- 
balance relaxation rate is much smaller than the elastic collision rate, and the conductivities 
(^) depend on the relation between frequency to and ^(A/T) 2 . 

In the frequency range u ^> ui(A/T) 2 in s-wave superconductors, where the factor 7 is 
real (see ]TT[), the weakly damped Carlson- Goldman mode appears. In c?-wave supercon- 



ductors the factor 7 = (7rAo/2T)^ii , i/uj) contains a large imaginary part due to the larger 
imbalance relaxation rate, and the related mode is highly damped. 

In the static limit our equations determine the penetration length Ie of the electric field 
into a ci-wave superconductor in direction a, when a current flows through a contact with a 
normal metal. Very near T c when one can neglect Andreev reflection of the quasiparticles 



we obtain l E = J (irA D a Ti)/(4:T), which agrees with the results of Choi [[H]. Here D a are 
diffusion coefficients related to the conductivities a^ a by the relation D a n 2 = 4iiaN a . The 
anisotropy of l\ is proportional to the conductivity anisotropy. 

If the order parameter is not of the pure (i-type symmetry, but is close to it: (A(0)) 2 <C 
(A(0) 2 ), then the results of our calculations are qualitatively the same. The main distinc- 
tions appear in the different energy and angle dependencies of the quasiparticle relaxation 
times. 

In conclusion, we calculated the linear response of layered <i-wave superconductors by 
means of the kinetic theory. We found the conductivities determining the quasiparticle 
currents created by the longitudinal and by the transverse electromagnetic fields. These 
results were applied to describe collective modes and the decay length of the electric field 
near the boundary with a normal metal. 

We are indebted to U. Eckern for reading the manuscript and helpful comments. 
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